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Micromechanics  and  Response  of  Particulate-Matrix  Fiber-Reinforced  Functionally 

Graded  Composites 


Abstract 

Reinforcement  of  fibrous  composites  by  stiff  particles  embedded  in  the  matrix  offers  the 
potential  for  simple,  economic  functional  grading,  enhanced  response  to  mechanical  loads  and 
improved  functioning  at  high  temperature.  The  solution  shown  in  the  paper  extends  a  version  of 
the  Mori-Tanaka  micromechanics  to  the  case  of  a  fiber-reinforced  material  with  inclusions 
embedded  within  the  matrix.  Furthermore,  we  present  several  bounds  and  estimates  for  the 
mechanical  properties  of  such  composite  and  show  how  functional  grading  within  these  bounds 
can  affect  the  performance  of  a  structure.  The  representative  example  presented  for  a  rectangular 
simply  supported  panel  subject  to  uniform  transverse  pressure  illustrates  an  abrupt  decrease  in 
deformations  due  to  the  addition  of  a  modest  amount  of  particles.  The  other  example  suggests 
superior  blast  resistance  of  the  panel  achieved  at  the  expense  of  only  a  small  increase  in  weight. 

1.  Introduction 

Hybrid  composite  materials  consisting  of  an  isotropic  polymer  matrix  reinforced  by  both 
particles  and  unidirectional  fibers  offer  the  potential  for  simple  functional  grading  to  tailor 
mechanical  response  and  reduce  stress  concentrations  around  attachments  and  discontinuities. 
Local  particle  reinforcement  can  increase  stiffness  and  strength  at  key  locations  at  the  expense  of 
a  relatively  small  increase  in  weight.  Moreover,  polymeric  composite  structures  subject  to 
thennal  loading  exhibit  matrix  deterioration  at  high  temperature,  but  a  matrix  reinforced  with 
ceramic  particles  offers  the  potential  to  increase  the  structural  endurance  and  load  capacity  at 
high  temperature  (Birman  and  Byrd,  2007). 
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The  application  of  hybrid  particulate-matrix  fiber-reinforced  composites  involves  an 
analysis  of  the  stiffness  of  such  materials  that  can  be  applied  to  structures  with  either  uniform  or 
variable  property  distribution.  For  example,  the  latter  structures  may  represent  thin-walled  shells 
or  plates  with  in-surface  variable  stiffness.  A  particularly  simple  and  attractive  grading  scheme 
involves  embedding  particles  only  in  the  outer  layers  of  a  laminate,  achieving  maximal  increases 
in  bending  stiffness  with  a  minimum  of  additional  weight.  Tailoring  of  the  volume  fraction  of 
particles,  especially  on  the  outer  laminae  of  a  composite,  is  far  simpler  than  the  alternative  of 
varying  fiber  volume  fraction  or  their  orientation. 

A  broad  range  of  homogenization  methods  exist  to  predict  the  properties  of  composite 
materials  including  particulate  inclusions.  Several  recent  reviews  and  articles  outlining  these 
methods  are  available  (e.g..  Tucker  and  Liang,  1999  and  Kakavas  and  Kontoni,  2006). 
Additionally,  the  study  of  Hu  and  Weng  (2000)  is  of  interest  as  it  outlines  and  compares 
micromechanical  models,  including  the  double-inclusion  method  (Hori  and  Nemat-Nasser  1993), 
and  the  models  of  Ponte  Castaneda  and  Willis  (1995)  and  of  Kuster-Toksoz  (1974)  and  the  well 
known  Mori-Tanaka  model  (1973). 

A  general  approach  to  the  characterization  of  a  hybrid  composite  consisting  of  three 
different  phases  was  proposed  by  Kanaun  and  Jeulin  (2001),  based  on  the  effective  field  method. 
This  method  is  based  on  the  assumption  that  the  strain  field  acting  on  every  inclusion  varies  for 
different  populations.  The  disadvantage  of  the  method  is  related  to  its  high  dependency  on  the 
precision  of  the  micromechanical  information  about  constituent  phases,  so  that  while  the 
theoretical  solution  may  be  accurate,  practical  implementations  would  depend  on  the  tight 
control  of  the  manufacturing  process  that  cannot  be  achieved  at  a  reasonable  cost  in  industry. 

The  analysis  presented  here  is  based  on  the  use  of  the  modification  of  the  Mori-Tanaka 
method  (Mori  and  Tanaka,  1973)  suggested  by  Benveniste  (1987).  In  addition,  the  upper  and 
lower  bounds  of  the  material  constants  (moduli  of  elasticity,  and  shear,  the  bulk  modulus  of  a 
particulate  material  and  the  Poisson  ratios)  are  determined  using  the  bounding  techniques. 
Available  bounds  include  those  by  Voight  and  Reuss  (Hill,  1952),  the  Hashin  method  for  a 
particulate  material  (Hashin,  1962)  combined  with  the  Hashin-Shtrikman  bounds  for  a  fiber- 
reinforced  material  (Hashin  and  Shtrikman,  1962,  1963),  the  Weng  approach  (Weng,  1992)  and 
the  three-point  bounds  technique  (Milton  and  Phan-Tsien,  1982).  The  accuracy  of  these  methods 
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has  been  debated  in  literature.  For  example,  the  modulus  of  elasticity  of  a  particulate  material 
can  be  accurately  evaluated  using  the  Mori-Tanaka  method  only  if  the  volume  fraction  of  the 
inclusions  does  not  exceed  40%  (Kwon  and  Dharan,  1995;  Sun  et  ah,  2007).  The  Hashin  bounds 
are  also  imperfect  predicting  an  exceedingly  wide  spectrum  of  the  modulus  of  elasticity  at  large 
particle  volume  fractions.  At  the  same  time,  the  predictions  obtained  by  the  Mori-Tanaka  and 
Hashin  methods  for  the  Poisson  ratio  of  particulate  composites  were  found  quite  accurate  (Sun  et 
ah,  2007).  On  the  other  hand,  Noor  and  Shah  (1993)  showed  that  the  Mori-Tanaka  method 
provides  an  accurate  prediction  of  the  properties  of  fiber-reinforced  composites  even  at  a  high 
volume  fraction  of  fibers.  In  the  following  analysis  the  volume  fraction  of  particles  is  assumed 
quite  small  reflecting  the  presence  of  fibers  that  occupy  a  significant  volume  within  the  material. 
Accordingly,  it  is  acceptable  to  rely  on  the  Mori-Tanaka  method  while  the  bounding  techniques 
are  also  applied  to  reflect  a  range  of  anticipated  variations  of  the  materials  properties. 

The  paper  presents  a  formulation  of  micromechanics  for  a  particulate-matrix  fiber- 
reinforced  material  based  on  the  extrapolation  of  the  Benveniste  approach  (1987)  for  the  case  of 
dissimilar  inclusions  (fibers  and  particles).  Subsequently,  the  bounding  techniques  outlined 
above  are  applied  in  a  two-step  solution:  the  bounds  are  first  established  for  a  particulate  matrix 
and  subsequently,  these  bounds  are  employed  to  establish  bounds  for  this  matrix  reinforced  by 
unidirectional  fibers.  Furthermore,  the  results  of  the  homogenation  conducted  by  on  the 
particulate-matrix  laminae  are  applied  to  predict  the  properties  of  a  cross-ply  material  consisting 
of  such  laminae.  Finally,  the  advantages  available  embedding  a  small  amount  of  particles  in 
fiber-reinforced  materials  are  illustrated  on  the  examples  of  a  simply  supported  cross-ply 
composite  plate  subject  to  static  pressure  or  dynamic  blast  overpressure.  In  both  these  examples, 
adding  glass  particles  to  a  glass/epoxy  plate  material  resulted  in  significant  reductions  of 
maximum  deformations  achieved  with  a  very  modest  additional  weight. 

2.  Benveniste-Type  Estimate  of  the  Stiffness  of  a  Fiber-Reinforced  Particulate  Material 

Consider  a  material  where  two  different  types  of  isotropic  inclusions  are  distributed 
within  an  isotropic  matrix.  The  properties  of  the  matrix  are  identified  in  the  subsequent  solution 
with  the  subscript  i  =  1,  while  two  types  of  the  inclusions  are  denoted  by  i  =  2  and  i  =  3.  We 
take  phase  2  to  be  spherical  particles  and  phase  3  to  be  aligned  fibers.  Each  type  of  inclusion 
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possesses  a  definite  geometry.  Then  the  stiffness  tensor  of  the  material  can  be  derived  as  a 
generalization  of  the  Mori-Tanaka  approach.  In  the  present  paper  this  generalization  is  based  on 
the  solution  by  Benveniste  (1987)  for  a  particulate  composite  with  a  single  type  of  inclusions. 
The  approach  is  based  on  the  following  assumptions: 

1.  All  material  phases  are  isotropic  and  linearly-elastic. 

2.  The  perturbed  strain  in  the  matrix  due  to  the  presence  of  inclusions  is  not  affected  by  the 
interaction  of  two  types  of  inclusions.  In  other  words,  each  type  of  inclusion  i  =  2  and 

i  =  3  affect  the  strains  in  the  matrix,  but  the  perturbed  matrix  strain  due  to  the  interaction 
between  these  inclusions  that  is  assumed  to  be  of  the  second  order  is  neglected. 

3.  Phase  3  is  represented  by  unidirectionally-oriented  fibers  of  circular  cross  section;  i.e., 
the  composite  material  is  a  lamina  with  embedded  particles.  This  assumption  is  only 
needed  to  utilize  the  Eshelby  tensor  for  cylindrical  inclusions.  In  general,  the  derivation 
shown  below  is  independent  of  the  orientation  and  shape  of  fibers  and  particles  as  long  as 
the  corresponding  Eshelby  tensor  is  known. 

The  average  stress  and  average  strain  tensors  for  the  material  under  consideration  are 
related  through  the  effective  stiffness  tensor: 

a-Ls  (1) 

where  the  bar  denotes  a  volumetric  average. 

The  effective  stiffness  tensor  for  a  matrix  with  two  different  types  of  embedded 
inclusions  can  be  given  as  a  generalization  of  the  expression  proposed  by  Hill  (1963): 

L  =  Li+^(L2-L1^A2+73^3_L1^A3  (2) 

where  L;  is  the  stiffness  tensor  of  the  /lh  phase,  and /j  are  volume  fractions  of  the 

corresponding  types  of  inclusions,  while  the  tensor  of  concentration  factors  A2  and  A3  represent 
the  relationship  between  the  tensors  of  average  strains  in  the  corresponding  inclusions  (A)  and 
the  mean  remote  strain  tensor  (soy. 

s2  =  A^q  ,  £3  =  (3) 
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Note  according  to  the  second  assumption,  this  approach  does  not  explicitly  account  for  the 
interaction  of  different  types  of  inclusions.  Accordingly,  it  is  applicable  only  if  at  least  one  type 
of  inclusion  has  a  relatively  small  volume  fraction. 

Expanding  the  Mori-Tanaka  (1973)  ideas,  the  tensors  of  average  strain  in  the  matrix  and 
in  inclusions  are  represented  as 

S2  =  £()  +£2+£3+  s'2  (4) 

£ 3  —  £q  +  £2  £$  £^ 

where  st  are  tensors  of  perturbations  superimposed  on  the  average  strain  in  the  matrix  as  a  result 
of  the  presence  of  the  corresponding  inclusions,  and  s\  are  tensors  of  average  perturbed  strain  in 
the  inclusions  relative  to  the  tensor  of  average  strain  in  the  matrix. 

The  tensors  of  average  stresses  in  the  inclusions  can  now  be  expressed  in  terms  of  the 
stiffness  of  the  matrix: 

^2  (*o  +  £1  +  S3  +  S2  )  =  In  0  +  S2  +  +  S2  ~  £1  ) 

Lj  -\-  £■,  +^3  +  s2 )  =  Lj  +  s2  +  £3  +  s2  —  s2  ^ 

where  s*  are  tensors  of  average  correlation  strain  in  the  corresponding  type  of  inclusions.  These 
tensors  are  related  to  the  tensors  of  perturbations  in  the  inclusions  defined  above  by 

£*2  =  S  2xs'2  si  =  S3  (6) 

where  S,  1  are  fourth-order  Eshelby’s  tensors.  These  tensors  are  presented  in  the  Appendix  1  for 
the  cases  of  spherical  inclusions  and  for  infinitely  long  cylindrical  inclusions  (fibers). 

It  is  evident  from  (4)  and  (5)  that  the  tensors  of  perturbation  strains  can  be  expressed  in 
terms  of  the  tensors  of  average  strain  in  the  corresponding  inclusions  as 

s;  =  s,l;'(l,-l,)5,  (1  =  2,3).  (7) 

As  also  directly  follows  from  (4), 


5 


ei=ei-el. 


(8) 


It  will  be  necessary  for  the  subsequent  transformation  to  express  the  tensors  of  average 
strain  in  each  type  of  inclusions  in  terms  of  the  tensor  of  the  average  strain  in  the  matrix,  i.e.  to 
find  the  coefficient  tensors  in  the  equations 

=  (9) 

This  is  easily  accomplished  using  (7)  and  (8): 

T^p  +  S^^.-L^J1  (10) 

where  I  is  a  fourth-order  unit  tensor. 

The  tensors  of  concentration  factors  are  now  determined  expressing  the  tensor  of  the 
applied  strain  that  also  represents  the  average  strain  in  an  equivalent  homogeneous  material  in 
tenns  of  strains  in  the  constituent  phases  through  the  rule  of  mixtures: 

£0  ~  f\£  1  fl£2  (11) 

Using  eqns.  (3)  and  (9)  in  (1 1)  yields 

A,-  =  %  (/J + f2%  +  /3T3J1  (12) 

Using  the  concentration  tensors  as  given  by  (12)  in  the  tensor  of  effective  stiffness  (2) 
yields  the  solution: 

3 

TJ=u+YJfi^i-uyi{fMf2%+h%f  (i3) 

i=2 

The  computations  utilizing  (13)  can  be  performed  using  the  tensor  decomposition 
presented  by  Walpole  (1983)  and  summarized  recently  by  Sevostianov  and  Kachanov  (2007). 
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Note  that  in  case  of  a  single  type  of  inclusions  this  result  converges  to  the  formula  derived  by 
Benveniste  (1987). 


3.  Two-Step  Approach  to  the  Evaluation  of  the  Elastic  Response  of  a  Particulate-Matrix 
Fiber-Reinforced  Composite  Material 

According  to  the  approach  adopted  in  this  paper,  conventional  bounds  and  estimates  are 
applied  to  evaluate  the  tensor  of  stiffness  of  the  isotropic  particle-reinforced  matrix. 
Subsequently,  this  tensor  of  stiffness  is  used  to  obtain  bounds  and  estimates  for  the  stiffness  of  a 
composite  comprised  of  isotropic  fibers  embedded  in  this  isotropic  particle-reinforced  matrix.  In 
this  section,  we  concentrate  on  three-point  bounds  since  as  shown  in  numerical  examples  they 
appear  provide  the  tightest  range  of  possible  values  of  material  constants. 

3.1  Bounds  and  Estimates  of  the  Mechanical  Response  of  the  Particle  Reinforced  Matrix 

Beran  and  Molyneux  (1966)  and  McCoy  (1970)  obtained  three-point  bounds  on  the 
effective  bulk  modulus  and  the  effective  shear  modulus  for  two-phase  composites.  These  contain 
two  calculated  parameters,  Q  and  rj2  (additionally,  )j\  =  \  -  )]i,  £\  =  \  -  Cf)  that  characterize  the 
shape  and  distribution  of  the  two  phases  (e.g.  Torquato,  1991)  and  must  be  evaluated 
numerically.  Milton  and  Phan-Thien  (1982)  further  improved  the  McCoy  (1970)  shear-modulus 
bounds. 

The  three-point  Milton-Phan-Thien  (1982)  bounds  on  the  effective  shear  modulus  Ge  for 
isotropic  two-phase  composites  are  given  by: 
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in  which  denotes  a  weighted  average  ((g)  =  G1^1  +  G2^2>  (G)  =  G^  +  G2£2 , 

(G)  =G1^1  +G2772),  and  a  tilde  represents  a  reverse-weighted  average,  e.g.  (Gj  =  G2</>i  +  G^2. 

Here,  (j)\  and  (fh  are  the  volume  fractions  of  epoxy  and  (relatively  stiff)  spherical  particles, 
respectively,  within  the  particulate-reinforced  matrix;  these  relate  to  f\  and  /2  as  =  /2  /(/)  +/2), 
with  <j)\  =  \-  tjh. 

The  simplified  form  (Milton,  1981)  of  the  three-point  Beran-Molyneux  (1966)  bounds  on 
the  bulk  modulus  Ke  of  isotropic  two-phase  composites  is  given  by: 
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where  d  =  3. 

Calculation  of  £?  and  ?/2  is  computationally  expensive,  and  values  have  been  reported  for 
only  a  limited  number  of  microstructures  at  this  time.  These  values  range  between  0. 15^2<4'2<^ 
and  Q.5(jh<ri2<(j>2  (Torquato,  1991).  For  randomly-spaced  spherical  particles,  £2-0.211  (fh  and 
772~O.48302  (Torquato,  2001).  The  dense  packing  limit  for  spherical  particles  is  that  fa  cannot 
exceed  approximately  0.63. 

Approximations  lying  between  these  bounds  are  available: 
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where  d=  3,  k2X  =  — K' 
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3.2  Bounds  and  Estimates  of  the  Mechanical  Response  of  the  Fiber-  and  Particle-Reinforced 
Lamina 

We  now  study  the  mechanics  of  a  lamina  containing  aligned  fibers  of  volume  fraction  fi 
embedded  in  an  isotropic  matrix  of  volume  fraction  f  A  f2  with  the  above  characteristics.  The 
dense  packing  limit  requires  f<  0.83.  Five  separate  moduli  are  needed  to  fully  describe  the 
linear  elastic  response  of  this  transversely  isotropic  material  but  only  three  of  the  effective 
moduli  are  independent  (Hill,  1963). 


Transverse  shear  modulus 

The  Silnutzer  three-point  lower  bound  on  the  transverse  shear  modulus,  Gj,  is  given  by 
Eq.  (14)  and  Eq.  (15)  (Silnutzer,  1972;  Torquato,  2001)  with  (j)x  replaced  with  (fx  +  fi),  (fa, 
replaced  with /3,  77 2  =  0.276 /3,  =  0.691  /3  +  0.0428  /32,  and  (  )  replaced  with  (  }  ,  where,  for 

example,  (G^  =  (/,  +  fffGpm  +  fGfib .  Here,  a  subscript  of fib  refers  to  properties  of  the  fibers 
(i=3),  a  subscript  of  pm  refers  to  properties  of  the  particle  reinforced  matrix  (i=l  and  2). 

The  Gibiansky-Torquato  (1995)  upper  bound  for  G],  is  tighter  than  the  Silnutzer  bound. 
Making  the  above  substitutions,  the  Gibiansky-Torquato  upper  bound  is  obtained  from  Eq.  (14) 
using  the  following  definition  for  0: 
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where  Kmax  is  the  greater  of  Kpm  and  Kflb,  Gfib>Gpm,  + 

//.^/l//2(67,;-67Jnj,and  Z  =  ^2(k£-K^  . 

An  estimate  for  Gj  within  these  bounds  is  provided  by  Eq.  (19),  with  d= 2. 


Transverse  bulk  modulus 


The  Silnutzer  three-point  lower  bound  on  the  transverse  bulk  modulus,  KTe,  is  given  by 
Eq.  (14)  and  Eq.  (15)  (Silnutzer,  1972;  Torquato,  2001)  with,  as  above,  (j)\  replaced  with  (/i  +  fi), 
(jh  replaced  with /3,  72  =  0.276  Qi  =  0.691  /3  +  0.0428  /32,  and  (  )  replaced  with  (  V,. 

A  tighter  upper  bound  is  the  Gibiansky-Torquato  upper  bound  (Gibiansky  and  Torquato, 
1995;  Torquato,  2001): 


P\  ,  GfGpm  +Kf 
'f  Kf  +  lG 


(21) 


An  estimate  for  the  lower  boundkij  is  obtained  by  Eq.  (19),  with  d= 2. 


Longitudinal-transverse  shear  modulus 

The  shear  modulus  for  distortion  of  the  laminate  in  axes  with  one  direction  aligned  with 
the  fiber  axes  is  given  by  the  Hashin-Rosen  (1964)  bounds: 
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The  Tsai-Halpin  semi-empirical  equations  provide  estimates  within  these  bounds  for  values  of 
the  parameter  ^  in  the  range  0  <  g<25: 
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The  value  E,=  2  has  been  shown  to  provide  a  good  estimate  for  laminae  containing  regularly 
spaced,  aligned  fibers,  and  will  be  adopted  in  the  following. 


Longitudinal  stiffness 

As  described  by  Torquato  (2001),  the  longitudinal  stiffness  Elc  of  a  transversely  isotropic 
two-phase  lamina  is  dictated  by  the  bounds  or  estimates  of  the  aforementioned  properties  (Hill, 
1964): 


E"  =  E)f +7 


(f pm  Vf) 


V  k  fib  kpm  J 


f  Kl  +  Gl/3, 


(24) 


where  vis  Poisson’s  ratio  of  the  isotropic  particular  matrix  (pm)  or  fibers  (f),  and  kt  =  K,  +  Gil  3 
for  each  phase  i. 

Lateral/transverse  Poisson ’s  ratio 

The  effective  Poisson’s  ratio  v[  T  is  dictated  by  the  other  material  constants  of  the 
laminate  (Hill,  1964;  Torquato,  2001): 
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(24) 


(V  pm  Vf) 


k 

V  pm 


KTe+GTJ  3y 


KfibJ 


where  vis  Poisson’s  ratio,  and  kj  =  Kt  +  Gj/  3  for  each  phase 

3.2  Moduli  and  bending  stiffness  of  a  0/90  laminate  composite 

Three  elastic  constants  are  required  to  assemble  the  effective  in-plane  mechanical 
properties  of  a  symmetric  0/90  laminate  with  equal  numbers  of  laminae  in  the  0°  and  90° 
directions.  These  can  be  written,  for  axes  parallel  to  the  fiber  directions,  as: 


_  2c12 

cn+c22  +2c23 


Eq  -  Cn+Cl2  -2cnv0 
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where  cn  =aEL 
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4.  Effectiveness  of  Particulate-Matrix  Fiber-Reinforced  Composites 

The  following  representative  examples  illustrate  the  advantages  available  adding  a 
relatively  small  amount  of  spherical  particles  to  the  matrix  of  fiber-reinforced  composites.  Two 
cases  considered  here  include  simply  supported  panels  subject  to  an  instantaneously  applied 
uniform  pressure  and  the  panels  undergoing  overpressure  due  to  an  explosive  blast.  We  will 
focus  only  on  choosing  the  volume  fraction  of  spheres  that  optimizes  weight  and  stiffness  of  the 
structure.  Cross-ply  panels  considered  here  have  dimensions  a  by  b  (in  the  x-  and  y-directions, 
respectively),  thickness  h,  and  density  p.  The  mass  m  of  such  a  panel  is 

m  =  pabt  (26) 
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The  peak  deflection  8  due  to  an  instantaneously  applied  uniform  pressure,  p,  on  one  face 
of  the  panel,  is  obtained  applying  a  magnification  factor  of  p  =  2  to  the  static  one  degree  of 
freedom  approximate  solution  that  was  shown  to  yield  an  error  equal  or  smaller  than  2.1%  in 
isotropic  panels  (Lekhnitskii,  1968): 
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plate  bending  stiffness  as  follows: 
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In  (28),  6  is  the  lamination  angle  of  the  layer,  0(/  are  reduced  stiffnesses,  Ei  and  E:  are  the  layer 

moduli  of  elasticity  in  the  fiber  direction  and  in  the  direction  perpendicular  to  fiber,  respectively, 
Gn  is  the  layer  in-plane  shear  modulus,  and  v12,  v2l  are  the  Poisson  ratios. 

In  addition  to  the  response  to  an  instantaneously  applied  pressure,  it  is  also  possible  to  use 

the  solution  of  Lekhnitskii  shown  above  to  predict  the  response  of  the  panel  to  an  explosive 
blast.  Problems  of  response  of  structures  to  blast  have  been  considered  by  a  number  of  authors, 
including  Houlston  et  al.  (1985),  Gupta  (1985),  Gupta  et  al.  (1987),  Binnan  and  Bert  (1987), 
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Librescu  and  Nosier  (1990),  Librescu  and  Na  (1998a,  1998b).  The  blast  overpressure  is  usually 
uniformly  distributed  over  the  surface  of  the  panel  and  can  be  presented  by  the  Friedlander 
exponential  decay  equation: 


p{t)  =  p  1-  — 


\ 

—  exp(-  at  I 1) 

tp) 


where  is  a  constant,  t  is  time,  t  is  a  positive  phase  duration  of  the  pulse  and  a' is  an 
experimental  decay  parameter. 

The  response  can  now  be  evaluated  using  a  single-degree  of  freedom  Lekhnitskii’s 
approximate  solution  combined  with  the  convolution  integral: 

S(t)  =  Sj(t)  (30) 

where  the  function  of  time  is  obtained  using  integrals  from  the  web  site  mathematica.com  as 
follows  (see  Appendix  2  for  a  justification  of  the  solution): 


/M=®J  1  -  exP  _  ~~  sink(<-S  =  t-»-T - [/(/)- /(())] 

o{  tp)  {  ip)  [(a)  +WP)  \ 

l{t)  =  exp  -y-  (o)Jp )[-  2  atp  +  ((a')2  -  (cantp  )2  \tp  - 1) 

/(°)  =  tp  {(0ntp  [-  2a'  +  (o')2  +  (contp  )2  ]cos  r oj  +  [(o')2  +  (contp  )2  \l  -  o')sin  a)nt\ 


-a  sin  (oj 


In  (3 1),  con  is  the  natural  frequency.  The  latter  frequency  could  be  detennined  using  (30)  in  the 

Ritz  solution  to  the  problem  of  free  vibrations  of  the  plate.  However,  considering  the  smallness 
of  the  error  obtained  using  the  static  Lekhnitskii  solution,  it  is  possible  to  use  the  exact  solution 
for  the  fundamental  frequency  of  the  plate,  so  that 

,  kUd .id.  A) 

m  \  a  a  b  b 


where  m  is  the  mass  of  the  plate  per  unit  surface  area. 
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5.  Numerical  Results 


Material  constants  of  an  epoxy  matrix  with  spherical  glass  inclusions  were  evaluated 
using  the  Mori-Tanaka  theory  and  compared  to  several  available  bounds  in  Fig.  1.  The  bounds 
included  in  this  figure  are  those  by  Voight  and  Reuss,  Hashin-Shtrikman,  and  three-point 
bounds.  Additionally,  the  weak  contrast  estimate  is  provided  in  Fig.  1.  Remarkably,  the  Mori- 
Tanaka  estimate  is  always  within  the  Voight-Reuss  bounds.  This  estimate  was  also  found  close 
to  the  lower  three-bound  estimate  by  the  Hashin-Shtrikman  and  by  three-point  bounds,  for  the 
elasticity,  shear  and  bulk  moduli,  while  the  prediction  for  the  Poisson  ratio  was  close  to  the  upper 
bounds  according  to  these  techniques.  It  should  be  noted  that  while  the  Mori-Tanaka  predictions 
remain  within  the  bounds  according  to  Hashin-Shtrikman,  they  were  slightly  outside  the  bounds 
by  the  three-point  technique,  particularly  at  larger  particle  volume  fractions,  though  the  deviation 
remained  quite  small.  Overall,  the  results  shown  in  Fig.  1  support  the  previously  cited  opinion 
that  the  Mori-Tanaka  theory  is  acceptable,  at  least  for  the  particle  volume  fraction  below  40%. 

The  results  shown  in  Fig.  2  refer  to  the  properties  of  a  particulate -matrix  fiber  reinforced 
material  (glass  fibers  and  particles  and  epoxy  matrix).  The  bounds  shown  in  this  figure  rely  on 
strict  three-point  bounds  from  Fig.  1  that  are  applied  to  the  case  of  a  fiber-reinforced  composite 
according  to  the  equations  shown  in  the  paper.  All  moduli  are  normalized  by  those  of  the  epoxy 
matrix.  As  follows  from  the  figure,  the  “Generalized  Benveniste”  estimate  corresponding  to  Eq. 
(13)  remains  within  the  bounds  only  if  both  the  fiber  as  well  as  the  particle  volume  fractions  are 
small.  It  is  interesting  to  note  that  while  the  present  theory  predicts  the  longitudinal  elastic 
modulus  that  is  higher  than  the  upper  three-point  bound,  the  transverse  modulus  is  lower  than  the 
lower  bound.  This  means  that  while  the  present  model  overestimates  the  longitudinal  stiffness  of 
a  particulate-matrix  fiber-reinforced  lamina,  the  transverse  stiffness  is  underestimated. 
Accordingly,  the  shortcomings  of  the  model  for  a  single  lamina  (that  is  never  contemplated  in 
realistic  design,  anyway),  may  cancel  each  other  in  a  cross-ply  material.  This  assumption  is 
further  investigated  in  Fig.  3.  Remarkably,  the  result  shown  in  this  figure  indicates  that  the 
predictions  for  the  stiffness  obtained  by  the  present  theory  for  a  cross-ply  material  remain  within 
the  tightest  three-point  bounds  (actually,  almost  coincide  with  the  lower  bound).  This  means  that 
the  theory  can  successfully  be  applied  to  the  prediction  of  the  stiffness  of  particulate -matrix 
fiber-reinforced  cross-ply  composites. 
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Appendix:  Eshelby’s  tensors  (Zhao  and  Weng,  1990) 

(1)  Components  of  Eshelby’s  tensor  for  a  spherical  inclusion  embedded  within  an  isotropic 
matrix: 


^1111  —  $2222  ~  ^3333 
^1 122  =  ^2233  =  ^3311 


*^1212  —  “^2323  —  ^3 131 


7-v, 

15(l-n) 

5vl  - 1 
15(1  -n) 
4  -  5vl 
15(l  -Vj) 


(Al) 


where  v,  is  the  matrix  Poisson  ratio. 


(2)  Eshelby’s  tensors  for  a  cylindrical  inclusion  embedded  within  an  isotropic  matrix 


Sun 

^2222 

c 

‘-’2211 

c 

lj2323 

C 

lj1212 


-9  -  V  -  0 

J1122  ‘-’1133  u 

_  C  _  5“4l/i 

_  3333  “8(1- v.) 

3311  2(1  -  Vj ) 

__  3  -  4vl 

“8(l-G 

=  S1313  =  0.25 


(A2) 


19 


Appendix  2:  Response  of  the  structure  modeled  by  a  single-degree-of- freedom  system  to  an 
arbitrary  dynamic  load  using  the  static  solution 

The  dynamic  response  of  a  single  degree  of  freedom  undamped  system  to  an  arbitrary 
load  is  governed  by  the  following  equation  of  motion  (all  notations  are  standard  for  the  theory  of 
vibrations): 

mw  +  kw  =  f{t)  (B 1) 

The  convolution  integral  describing  steady  state  vibrations  is 


w  = 


1 

mcon 


J  /(“=) sin  a>n{t  -  % )cU; 
0 


(B2) 


where  con  is  the  natural  frequency. 

Consider  a  structural  system  modeled  by  a  single-degree-of-freedom  approximation  as 
was  done  in  the  solution  of  Lekhnitskii  for  a  rectangular  simply  supported  plate.  Let  the  equation 
of  static  equilibrium  be 

L(S,)=p  (B3) 


p  being  a  static  load  and  Sst  the  static  solution  that  satisfies  (B3). 

If  the  same  system  is  subject  to  a  dynamic  load  pf{t),  the  equation  of  motion  is 
L(S)  =  pf(t)-mS  (B4) 

Assume  the  solution  in  the  form  8  -  8sIF{i  )  . 

Then  the  substitution  of  (B3)  into  (B4)  yields 
mSstF  +  L{Sst)F  =  fL{Sst )  (B5) 


F  = 


From  the  comparison  with  (Bl)  and  (B2)  it  is  obvious  that 

ran 


mS, 


'  J /fe) sin  o)n(t  ~  g)dg  =  ®„|/(^)sin  cojt  -  g)dg 


(B6) 


si  0 


where  the  natural  frequency  is  coH  = 


(ran 

mS,t 
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Figure  1.  Estimates  and  bounds  for  the  stiffening  of  an  isotropic  epoxy  matrix 
by  spherical  glass  inclusions. 
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Figure  2.  Estimates  and  bounds  for  the  stiffening  of  a  composite  containing  varying 
volume  fractions  /3  of  isotropic  glass  fibers  by  a  volume  fraction^  of  spherical  glass  inclusions. 
The  upper  and  lower  bounds  pictured  are  calculated  using  the  upper  and  lower  bounds  from 
Figure  1  in  conjunction  with  the  bounding  procedures  described  in  the  text.  Moduli  are 

normalized  by  those  of  epoxy. 
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Figure  3.  Estimates  and  bounds  for  the  stiffening  of  a  symmetric  0/90  cross-ply 
composite  containing  varying  volume  fractions  f$  of  isotropic  glass  fibers  by  a  volume  fraction^ 
of  spherical  glass  inclusions.  Moduli  are  normalized  by  the  elastic  modulus  of  epoxy. 
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